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^ Abstract 



The obtention (up to five or six loop orders) of nonasymptotic critical behavior, 
above and below Tc, from the field theoretical framework is presented and discussed. 

When talking about critical behavior one usually thinks of critical exponents (power 
laws), and eventually of corrections to scaling, all notions strictly related to the unprecise 
definition of an asymptotic critical domain. In fact criticality may be observed beyond 
that theoretical domain and, sometimes, this makes it difficult to compare theory and 
experiments. [Q] For example, it is thought that some systems could undergo a retarded 
^ ■ crossover from classical to Ising-like critical behaviors. In such a case, the critical 
{Sj i domain would be much larger than for, say, pure fluids. Consequently many correction- 
^ ! to-scaling terms should be introduced and, it is very likely that the series would not 
converge. For that reason, nonasymptotic theoretical expressions of critical behaviors are 
**^ ' required to describe such systems. 

*f • It is not very well known that, beyond the estimations of the critical exponents, the 
^h. renormalization group (RG) theory || is also adapted to provide us with nonasymp- 
r-| ! totic forms of the critical behavior especially when a crossover phenomenon occurs (the 
j>! \ crossover is then characterized by the competition of two fixed points). 

We briefly present here the principles of the calculations done within the massive field 
theoretical framework in three dimensions (d — 3) M and which have yielded accurate 
nonasymptotic forms of the susceptibility x{ T ) an d the specific heat C(r) for r = (T — 
T c )/T c > and r < 0, of the correlation length £(r) for r > and of the coexistence curve 
M(r) for r < 0. [|], H The calculations presented here have induced, directly or indirectly, 
several subsequent works. ]7j || We hope that this text will encourage further works on 
nonasymptotic critical behavior. We think, in particular, that the variational perturbation 
theory used recently to estimate universal exponents |J and amplitude ratios, |IJ could 
be an advantageous tool. 

Let us first specify the meaning of the title. "Nonasymptotic critical behavior" means 
that we perform a resummation of the infinite series of correction-to-scaling terms which 
are expected [JTTJ] in the asymptotic expression of any singular quantity such as £(t). 
Particularly, for r — > + ~, we have: 



£{t) = \r\ 
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in which v is a critical exponent, £g stands for the leading critical amplitudes in the two 
phases and the coefficients a£+,- correspond to the amplitudes of the confluent corrections 
to scaling controlled by the exponents A m (m = 1, 2, ■ • ■ , oo). Those exponents [y and 
A m ) arise in a linear study [II] of the RG transformation in the vicinity of a fixed point: 
the solutions of the eigenvalue problem provide us with some positive (say one, A , for 
simplicity^ and infinitely many negative eigenvalues (A m for m = 1, 2, • • • , oo). Then we 
have v = 1/Ao, A m = — u\ m (m = 1,2, •••,oo). The case m = 1 corresponds to the 
first correction-to-scaling term associated to the largest negative eigenvalue and the usual 
notations are ui (for — Ai) and A = uju (for Ai). In a linear study of the RG, discarding 
the next-to- leading correction terms in Eq. (p]), it is usual to write: 



i{r) ~ \t\ 



1 + a 



T 



(2) 



with the universality of the ratios £o~/£o~ and a^/aj. Eq. (|2]) is only valid asymptotically 
close to the critical point. 

Notice that the infinite sum in Eq. ([[]) does not converge for large values of r. Thus, 
to get a useful nonasymptotic expression of the critical part of £, we must consider a 
resummation procedure. It is provided to us by the RG theory. However the framework 
we use implies an approximation: "from field theory" means that only one family of 
correction-to-scaling terms (associated to m — 1) is accounted for. Thus, instead of Eq. 
(P, our effective expression for £(r) is: 



£ft(t) = & |t| 



1 + 



.(») 
n=l 



in A 



(3) 



As one could see from the use of our calculations in analyzing experimental [O] or Monte- 
Carlo ][14| data, the approximation of field theory does not prevent the study from yielding 
physically useful nonasymptotic critical behaviors. In fact, in the case of Eq. (|^), the range 
< r < oo corresponds to an interpolation between two fixed points and, consequently, 
the crossover is described by universal functions. 
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1 Principles of the calculations 



One starts from the "bare" or unrenormalized 4 -hamiltonian in d euclidean dimensions 
corresponding to the scalar field theory to be renormalized: 



H 



d d x 



{(V0 o ) 2 + r o ^} + !^ 



(4) 



in which it is implicit that 0o depends on x and may eventually represents a vector with 
n components in which case Eq. (|J) is supposed to satisfy the 0(n)-symmetry. The bare 
coupling g is dimensionful and "measured" in unit of A (the ultra-violet cutoff): 

5-o = u A £ 



lr The number of positive eigenvalues depends on the fixed point considered. Since we are interested 
in a critical point (we mean not a multi-critical point) then there is only one positive eigenvalue and the 
fixed point referred to in that case is the famous Wilson-Fisher p2| fixed point. 
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in which e = 4 — d and Uq is dimensionless. 

There exists two kinds of renormalization scheme for the scalar field theory: the mas- 
sive and the Weinberg |16|] schemes. In the massive scheme, the unit of reference is 
provided by the mass parameter m = In this framework the critical theory, corre- 
sponding to m — 0, is not defined. On the contrary, in the Weinberg scheme, one first 
defines the critical theory (massless theory) and the unit of length scale is provided by 
the inverse of some arbitrary momentum-subtraction-point parameter \x. In that renor- 
malization scheme, the "soft-mass" parameter t is introduced via the renormalization of 
insertions of the ^-operator within the vertex-functions; when it is different from zero, t 
is (linearly) proportional to the reduced temperature scale r defined above. Though we 
have effectuated our calculations within the massive framework — because the longest 
available^] perturbative series [[17] have been obtained within the massive scheme directly 
in d = 3 — , the presentation of the principles of the calculation is simpler within the Wein- 
berg scheme. In that scheme, the renormalization conditions correspond to the following 
(reVdefinitions: 



Z 3 (u 



u A e 



r 



/i u 



r 0c + 



Zi (u) 

[Zs («)] 2 

Z2(U) 

Z 3 (u) 



(5) 
(6) 

(7) 



in which u is the renormalized coupling and tq c is defined by: 



ro°' 2) (p; r ,go) 







(9) 



p=0,ro=ro, 

in which the subscript refers to the bare theory. 

Up to analytical terms which are usually neglected when studying critical phenomena 
the quantity tq — tq c is proportional to the physical parameter r: 



A 2 



0T + O r 



(10) 



in which 9 is a nonuniversal factor. 

The renormalized iV-point vertex-functions with L- insertions of the (0) 2 operator are 
related to their unrenormalized counter parts as follows: 



r( L ' N H{q,p};u^) = [z 3 (u)] N/2 



Z 2 {u) 
Z 3 (u) 



r 



(L,N) 



({q,p};r ,g ) 



That the renormalization functions (u) be defined by renormalization conditions 
on the 2-point and 4-point vertex-functions considered at some subtraction momentum 
point expressed in terms of the only dimensionful (momentum-like) parameter fi or by a 
"minimal" subtraction procedure (such as the subtraction of poles located at e = 0) does 
not matter for the following. 



2 T he unpublished Guclph report | 17] may be obtained via H. Kleinert and V. Schulte-Frohlinde web 
site at http: / /www. physik.fu-berlin.de/' kleinert /kleiner jeb8/programs/programs. html. 
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In field theory, the RG originates from the arbitrariness of the subtraction procedure 
for a given bare theory. Hence, the renormalized quantities u and t become functions of 
the renormalization parameter I = — In (///A) • Consequently, Eq. ([8]) must be understood 
as follows: 

r° = r 0e + pP§ m (11) 

Now, by imposing that t(l) remains a fixed quantity (say t (I) = 1), one relates the 
evolution of u (/) to the approach to the critical point of the bare (physical) theory (defined 
by r - > r 0c)- Then, with t (/) = 1, Eq. ([□]) shows that, for r = r 0c , u (/) must take on a 
particular value u* so that Z 2 (u*) /Z 3 (u*) vanishes. Of course, u* is the nontrivial zero 
of the famous /5-function: 

nu) = - du{l) 



di 



(12) 



with a positive u> = d/3 (u) /du\ u=u , so that u (I) — > u*. 

The pure scaling (power law) regime of vertex-functions corresponds to u (I) = u* and 
the first correction-to-scaling term [as in Eq. @] to a linear correction proportional to 
u (/) — it*. As u (I) moves further away from u*, more and more correction terms must be 
included but then a nonlinear study is required. It is a matter of fact that the domain 
< u (I) < u* corresponds to the entire domain oo > r — r 0c > 0. Therefore, if one 
re-sums perturbative series in powers of u (I) in the range ]0, u*[, one implicitly obtains 
nonasymptotic critical answers which interpolate Q between a classical critical behavior 
(when u (I) is small) and, say, an Ising-like critical behavior (for the O(l)-symmetry) 
when it (/) approaches u* . It remains to invert Eq. (|TI|) to express these answers under 
the forms of functions of r — r 0c [or of r, via Eq. ([K])] which is the genuine physical 
"measure" of the distance to the critical point. 

With a view to get the best possible accuracy, we have looked at the available cal- 
culations up to relatively high orders of perturbation series. There are two kinds of such 
calculations: 

1. analytically up to fifth order in the Weinberg scheme with dimensional regularization 
and minimal subtractions. fEDl 



2. numerically up to sixth order for d — 3 in the massive scheme. 17 



In the two cases, only the renormalization functions Z{ are considered. This is because 
the theoretical interest is usually focused on the critical exponents, the obtention of which 
comes from the knowledge of the Z^s. For example the series expansion for the critical 
exponent 77 are given by: 

v(u) = P(u)-^-\nZ 3 {u) (13) 

once considered at u = u*. 

However, we are not simply interested in the critical exponents but in complete func- 
tions such as £ and x- Now, only in the case 2, the renormalization conditions are such 
that £ and x are known in terms of the ZiS. This is not true in the case lfl. 



3 We have also performed calculations |l8| for u > u* . In this range, the sign of the first correction- 
to-scaling term is changed and corresponds to the Ising model . 

4 This is why the amplitude functions are known only up to three loop order in this scheme. @ 



4 



We denote the renormalized parameters of the massive scheme by g and m (instead of 
u and t). Their relations to the bare parameters are similar to those given by Eqs. (0-0) 
except that in addition to the change u — > g, Eqs. (|7|-[|) now read: 



9o 



r 



m g 



Zi (g) 

[Zs (g)f 



5m 2 + 



m 



z 3 (g) 



(14) 
(15) 



The mass shift Sm 2 is defined by a subtraction condition^ which avoids the explicit 
consideration of r 0c via Eq. (0), namely: 



2 -r(°' 2 ) (p;m,g) 



p=0 



r(°' 2 ) (0;m,g) = m 2 
The other subtraction conditions^] which define the Z^s read: 

_d_ 

dp 2 

T^({0};m,g) 
r (1 ' 2) ({0,0};m )5 ) 

so that the physical (bare) quantities £ and x are given by: 

[Zs (g)} 2 



(16) 



m g 
1 



r\g) 

x-\g) 



m = g - 



ry-l 2 

Z 3 m 



gZi (g) 



go 



[Zs (g)Y 



[gZi (g)] 2 

The re-summations of the perturbative series for those quantities have been done using 



the technique initiated by Le Guillou and Zinn- Justin HjJ after having accounted for the 
singularities of the Zj's at the fixed point g*. They may be easily treated by writing, e.g. 
for Z 3 (g) which has a singularity at g* of the form (g* — g) v . 



Z 3 {g) = Z 3 (y) exp 



in which y is some small value of g, the definitions of f3 (x) and rj (x) being unchanged 
in their forms compared to Eqs. 




13]). Let us mention that some difficulties could be 



encountered in the resummation procedure due to nonanalytic confluent singularities |22 
in the /5-function at g*, but to date they have not been numerically observed. 

Thus, in the homogeneous phase, the physical quantities £ and x can be easily esti- 



mated as functions of g in the range ]0,g*[ from the calculated series. JT^ However our 
aim was to obtain those quantities as function of r (i.e. of r — r 0c ). Now, the massive 
framework uses mocr" instead of t oc r, therefore the linear relation to r is lost. In order 
to reintroduce it, we use the fact that, at zero-external- momenta: 



5 Which eliminates the quadratic ultra-violet divergences occuring at d = 4. 
6 Which eliminate the logarithmic ultra-violet divergences occuring at d = 4. 
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[Z2 (g)} 



1 = T { q' 2) (r Q ,g ) 



_d_ 

Ore 



ri°' 2) {r ,g ] 



Using Eqs (H-pf), we reexpress this under the following form: 

d{r Q /gl) 



dg 



W dg \ [gZ, (g)f 



which, after integration allows us to (implicitly) define an effective (and nonperturbative) 
critical value r' 0c by reference to the fixed point value g*\ 



t(g) 



Oc 



gl 



dx Z2 (x) 



d j [Z 3 (x)f 



dx { [ X Z 1 (x)Y 



(17) 



The integrand of Eq. (|T7|) may be estimated using the same procedure as before and 
the integration has been done numerically yielding the numerical evolution of t (g) in the 
range ]0, g*[. The final results (the functions £(t) and x{t)) were obtained after a fitting 
procedure of the implicit form £(<?), x{g) an d t (g). This summarizes the calculations done 
in the homogeneous phase || which included also the specific heat C(t) the perturbative 
series of which were previously |23| extracted from the Guelph report; [|lj]] the calculations 
were performed for the symmetries n = 1, 2 and 3. 



2 Calculations in the inhomogeneous phase and the 
critical bare mass 

"...it is more difficult to calculate physical quantities in the ordered phase because the 
theory is parametrized in terms of the disordered phase correlation length m _1 which is 
singular at T c . Also the normalization of correlation functions is singular at T c ." |24] 



We have calculated || the perturbative series for the free energy directly at d = 3 
using the numerous already-estimated Feynman integrals of the massive scheme [17] and 
new kinds of integrals which have been estimated for the occasion. 

Because the free energy is generally written in terms of T — T c , we have been led to 
explicitly consider the delicate question of the critical bare mass. Indeed, it is known 
that the perturbative series of super-renormalizable massless field theories (such as <^ <4 ) 
develop infrared singularities which are usually simply ignored within the e-expansion 
framework. In 1973, using a dimensional regularization, Symanzik p5f has shown that 



the critical bare mass — which has the form ro c = g^ e f (e) in which / (e) has poles at 
e — 2/k (k — 1,2, ...00) — is in fact an infrared regulator for the theory. However the final 
result (free of infrared divergences) is no longer perturbative (e.g., logarithms of go appear 
at d — 3). Though the nonperturbative nature of r 0c is an important aspect of the RG 
theory, 0] this question may be circumvented when looking at the critical behavior. This 
is because T c is a nonuniversal quantity, thus its explicit determination is not required, 
only the difference T — T c is needed. Hence, provided that Eq. (||) be again satisfied, 
one may redefine r 0c [as it has been done in Eq. (PI)]. Consequently, it is allowed to 
perform a particular mass-shift r = r' Q + 5r (e) in such a way that 5r (e) subtracts the 
poles occuring at e = 1, and to fix afterwards the critical temperature in terms of r' . 
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The series for the free energy have then been obtained graph by graph up to five loops 
according to the following rules: 



1. graphs involving only 4 - vertices which were already estimated [T7j with the mass- 



shift parameter Sm 2 [defined by Eq. (|T^, 16)] have been re-expressed to account for 
the mass-shift parameter 5r (e) . 

2. New Feynman integrals at d = 3 with their weights involving: 

• exclusively </> 3 -vertices have been calculated and compared to existing esti- 
mates. [^7] 



• 3 -vertices mixed with a single 4 -vertex have been estimated for the first time 
for the occasion. 

Those series for the free energy have been used by Guida and Zinn- Justin [§] to give an 
accurate estimation of the scaled equation of state. But this kind of consideration does 
not account for any correction-to-scaling terms and the comparison with experiments 
is not easy. Instead, we were again interested in actually measurable quantities like 
the susceptibility \i the specific heat C and the spontaneous magnetization M in the 
inhomogeneous phase. We have not considered the correlation length £ in this phase 
because the required Feynman integrals were not calculated by Nickel et al; [17]] this 
quantity has been considered afterwards at d — 3 but up to 3-loop order only. fl28 | 




Because the renormalization procedure is unchanged in going into the broken-symmetry 
phase, the critical singularities at the fixed point g* may be taken into account with the 
same renormalization functions Zi (g) as defined previously. But the relation between 
the linear measure of the distance to T c and the unchanged renormalized coupling g is 
different. Instead of Eq. ( |iT|) we have obtained: 

t' (g) = - 

in which U (g) is given || as power series in g. 

Obviously, our nonasymptotic study of the critical behavior accounts for all the univer- 
sal properties which are expected when r — > 0. Consequently as a by-product, estimates 
of universal combinations of leading critical amplitudes were given for the first time from 
the five loop order at d = 3. The recent careful re-estimations || |29|] of those universal 
combinations from the same series have mainly reduced the error-bars. We also gave for 
the first time accurate estimates of some |3tJ of the universal ratios of the first confluent 
correction-to-scaling 0. 
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